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\ Abstract. The correlation between black hole mass Mbh and the cen- 

■ tral stellar velocity dispersion a in nearby elliptical galaxies affords a 
I \ novel way to determine Mbh in active galaxies. We report on measure- 

■ ments of a from optical spectra of 7 BL Lac host galaxies. The range of 
I \ a (160 - 290 km s~^) corresponds to estimated Mbh of 5 x 10'' - 9 x 10*^ 

■ Mq. The average ratio of Mbh to the host galaxy mass is 1.4 x 10~^, 
\ consistent with that estimated in other active and inactive galaxies. The 

■ range of a and Mbh of the BL Lacs are similar to those obtained for 
low redshift radio galaxies, in good agreement with the predictions of the 

O ' unified models for radio-loud active galaxies. 

■ 

^> : 

\ 1. Introduction 

H ; 

_Cr ! The mass of the central black hole (BH) is of paramount importance in theoret- 

ical models of AGN. In particular, the dependence of BH mass {Mbh) on global 
host galaxy properties provides clues to the role of BHs in galaxy formation 
and evolution. Dynamical determination of Mbh in AGN is difficult because of 
the bright nuclear emission. The main method that has proved to be successful 
for AGN is the time-consuming reverberation mapping of broad emission lines. 
Only for a few quasars and Seyfert galaxies Mbh is thus known (e.g. Kaspi et al. 
2000; Wandel 2002). Reverberation mapping cannot obviously be employed for 
BL Lac objects because they lack prominent broad emission lines. The discov- 
ery of a correlation between Mbh and the bulge luminosity in nearby early-type 
galaxies (e.g. Magorrian et al. 1998) offered a new tool for evaluating Mbh (see 
Merritt & Ferrarese 2001). This correlation has been applied so far for nearby 
quasars (McLure & Dunlop 2001) and BL Lacs (Treves et al. 2002). 

Recently, a tighter correlation was found relating Mbh with the stellar 
velocity dispersion a of the bulge in nearby inactive galaxies (Gebhardt et al. 
2000; Ferrarese &: Merritt 2000). This relationship clearly demonstrates a con- 
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nection between BHs and bulges of galaxies and has spurred substantial effort 
in theoretical modelling (e.g. Silk &: Rees 1998; Hahnelt & Kauffmann 2000). 
The relationship predicts more accurately Mbh, but requires the measurement 
of a in AGN host galaxies that is difficult, in particular for objects at high red- 
shift and with luminous nuclei. On the other hand, BL Lacs have relatively 
fainter nuclei than quasars, and for low redshift BL Lacs this measurement can 
be secured with a single spectrum observable with a medium-sized telescope. 

We have carried out medium resolution optical spectroscopy of the host 
galaxies of nearby (z < 0.2) BL Lac objects to derive a and to determine Mbh- 
Here wc present results based on a sample of 7 BL Lacs including 2 LBLs (3C 
371 and PKS 2201+04) and 5 HBLs (see also Falomo, these proceedings; and full 
discussion in Falomo, Kotilainen &: Treves 2002), For all of them, high quality 
images have been obtained either from the ground (Falomo k. Kotilainen 1999) 
or with HST (Urry et al. 2000; Falomo et al. 2000). From these images, the 
characterization of the host galaxies and the nuclear luminosity can be obtained. 
In Falomo, Kotilainen k, Treves (in prep.), we shall discuss the implications of 
the results of the full sample of 11 BL Lacs for the unified models of radio-loud 
AGN. 

2. Observations and data analysis 

The observations were obtained using the 2.5m Nordic Optical Telescope (NOT) 
equipped with ALFOSC. Spectra were secured using two grisms to cover the 
spectral ranges 4800 - 5800 A (setup A) and 5700 - 8000 A (setup B) at 0.54 
A pixel""^ and 1.3 A pixel""^ dispersion, respectively. This allows us to measure 
the absorption lines of H/3 (4861 A), Mg I (5175 A), Ca E-band (5269 A), Na I 
(5892 A) and the TiO + Cal (6178 A), TiO + Fel (6266 A) and other absorption 
line blends from the host galaxies at a spectral resolution R ~3000. 

The chosen grisms combined with a 1" slit yield a spectral resolution for a 
measurement of ~60 - 80 km s~^, which is adequate for the expected range of cr 
in luminous ellipticals (e.g. Djorgovski & Davis 1987; Bender, Burstein & Faber 
1992) such as the hosts of BL Lacs. In addition, we acquired spectra of bright 
stars of type G8-III to Kl-III, that exhibit low rotational velocity (Vxsin (i) < 
20 kms~^) to be used as templates of zero velocity dispersion. Furthermore, 
spectra of the well studied nearby elliptical galaxy NGC 5831 were secured to 
provide a test of the adopted procedure to derive a. 

During the observations, seeing ranged between 1" and 1.5". The slit was 
centered on-target or positioned 1" away from the nucleus and the ID spectrum 
was extracted from an aperture of 3" - 5" diameter, in all cases within the 
effective radius of the host galaxy. In one case (Mrk 180), spectra with the slit 
both on-target and off-centered by 1" were taken but no significant difference 
was found in the shape of the spectral features. 

The stellar velocity dispersion a was determined using the Fourier Quotient 
method (e.g. Sargent et al. 1977). The spectra were normalized by subtracting 
the continuum, converted to a logarithmic scale and multiplied by a cosine bell 
function that apodizes 10% of the pixels at each end of the spectrum. All 
regions affected by emission lines (sec Fig. 1) have been masked in the analysis. 
Finally, the Fourier Transform of the galaxy spectra were divided by the Fourier 



APS Conf. Ser. Style 



3 



Transform of template stars and a was computed from a fit with a Gaussian 
broadening function (see e.g. Bertola et al. 1984; Kuijken Sz Merrifield 1993). 
The rms scatter of the a results using different template stars was typically ~10 
kms~^ and can be considered as the minimum uncertainty of the measurement. 

For three objects we have spectra in both spectral ranges. The resulting 
values of a are in good agreement, ensuring sufficient homogeneity of data taken 
with different grisms and/or resolution. Note, however, that there is a tendency 
for the red grism (lower resolution) data to result in slightly larger value of a. For 
the nearby elliptical NGC 5831 we obtained a = 167=b5 and 185ibl0 km s~^ for 
the setup A and B, respectively, in good agreement with previous measurements 
(< o" > = 168 km s^^; Prugniel et al. 1998). In Fig. 1 we show the spectrum of 
the BL Lac object PKS 2201+04, observed in both spectral ranges. 




Figure 1. Optical spectrum of the BL Lac object PKS 2201+04 (z 
= 0.027), normalized to the continuum and plotted in the rest frame. 

Since early-type galaxies exhibit gradients in velocity dispersion (e.g. Fisher, 
Illingworth Sz Pranx 1995), the measured value of a depends somewhat on the 
distance of the galaxies and the aperture size (see Falomo et al. 2002 for details). 
The individual measurements of corrected a are given in Table 1. 

Barth, Ho & Sargent (2002) have recently reported a similar systematic 
study of (7 in 11 BL Lac objects, six of which are common with our sample. 
There are some discrepancies for individual BL Lacs, most notably for Mrk 
501 and I Zw 187, for which they derive a = 372+18 and 171+12 km s^S 
respectively. These differences are most likely related to the method of deriving 
a (Fourier Quotient vs. direct fitting) and the different wavelength range used. 
The average values of a, however, are in good agreement in the two samples. 



4 



Author <fe Co-author 



Table 1. Velocity dispersion and BH masses of BL Lacs. 



Object 


z 




log{MBH)a 


log{MBH)bulge 


log(M(, 






km s^-*^ 


[Mo] 


[Mo] 


[Mo] 


Mrk 421 


0.031 


236+10 


8.50+0.18 


8.65 


11.20 


Mrk 180 


0.045 


244+10 


8.57+0.19 


8.50 


11.45 


Mrk 501 


0.034 


291+13 


8.93+0.21 


9.00 


11.59 


I Zw 187 


0.055 


253+15 


8.65+0.18 


8.20 


11.39 


3C 371 


0.051 


284+18 


8.88+0.20 


8.90 


11.32 


lES 1959+65 


0.048 


195+15 


8.12+0.13 


8.30 


11.27 


PKS 2201+04 


0.027 


160+7 


7.72+0.13 


8.27 


11.00 



3. Results and discussion 



We have adopted the relationship between Mbh and o"c found for nearby early- 
type galaxies based on optical spectroscopy (Merritt &: Ferrarese 2001): Mbh = 
1.48+0.24 X 10^ (0-/200 km s-i)4-65±o.48 jjyt^] assume that this relationship 
is also valid for BL Lacs. This is consistent with imaging studies of BL Lacs 
(Falomo & Kotilainen 1999; Urry et al. 2000; Falomo et al 2000), indicating 
that all BL Lacs are hosted by giant ellipticals. The derived values of Mbh 
are reported in Table 1, where the errors are the composition in quadrature of 
uncertainties in a and in the Merritt &; Ferrarese (2001) relationship. Using 
instead the Gebhardt et al. (2000) relationship tends to yield slightly lower 
values of Mbh but does not substantially modify our main conclusions. The 
values of Mbh (Table 1) span a factor ~20 from 5 x lO'^ M© for PKS 2201+04 
to 9 X 10^ Mo for Mrk 501, with a median value of 4 x 10^ Mq. These values 
are in good agreement with those derived by Barth ct al. (2002). 

As mentioned above, Mbh is also correlated with the bulge luminosity of the 
host galaxy. Mbh was calculated following the relationship by McLure Sz Dunlop 
(2002): log Mbh = -0.50+0.05 Mr ~ 2.91+1.23 [Mq]. The corresponding values 
of Mbh are given in Table 1. For most sources the difference of Mbh derived 
with the two methods is within the estimated uncertainty. The average values 
of Mbh derived, respectively, from a and the host luminosity are: <logMBij >o. 
= 8.62 + 0.23 and <logMBH >host = 8.66 + 0.25. 

The measurements of a combined with the effective radii of the host galaxies 
can be used to estimate the mass of the hosts through the relationship (Bender 
et al. 1992): Mw=5cJ^re/G. This dynamical mass (Table 1) is in the range of 
1 - 4 X 10^^ Mq. The ratio between Mbh and M^ost is iii the range of 0.5 - 3.6 
X 10~^, with average < Mbh / M^ost > = 1-4 x 10~^. This is in good agreement 
with values derived for both AGN and inactive galaxies (< MBn/Mfiost > = 1-2 
x 10"^ ; McLure & Dunlop 2001; Merritt & Ferrarese 2001). 

According to the unified model of radio-loud AGN (e.g. Urry & Padovani 
1995), BL Lacs are drawn from the population of FR I radio galaxies viewed 
along the jet axis. It is therefore interesting to compare orientation-independent 
properties of BL Lacs and radio galaxies. The largest comparison sample with 
available measurements of a, effective radii and surface brightness is that of 73 
low redshift radio galaxies by Bettoni et al. (2001, 2002). In Fig. 2, we compare 
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the distribution of the BH masses of the BL Lacs and the radio galaxies from 
Bettoni et al. 
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Figure 2. The distribution of the Mbh of the BL Lac objects com- 
pared with that of low redshift radio galaxies (from Bettoni et al. 2002). 

The close similarity of the distributions in these samples implies that the 
distributions of Mbh in BL Lacs and radio galaxies are indistinguishable, consis- 
tent with both types of AGN belonging to the same parent population, observed 
from different orientation angles (see also Barth ct al. 2002). Additional support 
to this view is given by the similarity of BL Lacs and radio galaxies on the funda- 
mental plane (see e.g. Falomo, these proceedings; Barth et al. 2002). However, 
due to the small number statistics of the BL Lacs, and the incompleteness of the 
samples, this result can not yet be considered conclusive. Furthermore, there is 
a clear dearth of BL Lacs with the largest Mbh (> 10^ M0), consistent with 
BL Lac hosts being biased toward less massive and less luminous host galaxies 
than radio galaxies (e.g. Urry et al. 2000). 



6 



Author Sz Co-author 



References 

Barth,A., Ho,L.C., Sargent,W.L.W., 2002, ApJ, in press ( |astro-ph/0209562| ) 

Bender,R., Burstein,D., Faber,S.M., 1992, ApJ 399, 462 

Bertola,F., Bettoni,D., Rusconi,L., Sedmak,G., 1984, AJ 89, 356 

Bettom,D., Falomo,R., Fasano,G., et al., 2001, A&A 380, 471 

Bettoni,D., Falomo,R., Fasano,G., Govoni,F., 2002, A&A, in press 

Djorgovski,S., Davis,M., 1987, ApJ 313, 59 

Falomo,R., Kotilainen,J., 1999, A&A 352, 85 

Falomo,R., Scarpa,R., Treves,A., Urry,C.M., 2000, ApJ 542, 731 

Falomo,R., Kotilainen,J.K., Treves, A., 2002, ApJ 569, L35 

Ferrarese,L., Merritt,D., 2000, ApJ 539, L9 

Fisher,D., Illingworth,G., Franx,M., 1995, ApJ 438, 539 

Gebhardt,K., Bender,R., Bower, G., et al., 2000, ApJ 539, L13 

Hahnelt,M.G., Kauffmann,G., 2000, MNRAS 318, L35 

Kaspi,S., Smith,P.S., Netzer,H., et al., 2000, ApJ 533, 631 

Kuijken,K., Merrifield,M.R., 1993, MNRAS 264, 712 

Magorrian,J., Tremaine,S., Richstone,D., et al., 1998, AJ 115, 2285 

McLure,R., Dunlop,J., 2001, MNRAS 327, 199 

McLure,R., Dunlop,J., 2002, MNRAS 331, 795 

Merritt,D., Ferrarese,L., 2001, MNRAS 320, L30 

Prugniel,Ph., Zasov,A., BusareIlo,G., Simien,F., 1998, A&AS 127, 117 
Sargent,W.L.W., Schechter,P.L., Boksenberg,A., Shortridge,K., 1977, ApJ 212, 
362 

Silk,J., Rees,M.J., 1998, A&A 331, LI 

Treves, A., Carangelo,N., Falomo,R., et al., 2002, Issues in Unification of AGN, 

eds. R.Maiolino et al., p. 303 (ASP) 
Urry,C.M., Padovani,P., 1995, PASP 107, 803 
Urry,C.M., Scarpa,R., 0'Dowd,M., et al., 2000, ApJ 532, 816 
Wandel,A., 2002, ApJ, 565, 762 



